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ABSTRACT 

Coherent Anti-Stokes Raman Spectroscopy (CARS) is a new 
type of Raman Spectroscopy. The phenomenon is associated with 
the nonlinear conversion of two laser beams in a medium into a 
third collimated beam at the anti-Stokes frequency of the 
medium. This analysis evaluates the performance capabilities 
of CARS under various experimental conditions for the molecules 
O2/N2/H2 and CO. Laser intensity fluctuations and shot noise 
are introduced as noise sources to predict signal-to-noise 
ratios (S/N) . The S/N is evaluated as a function of the par- 
tial pressure of a gas to measure the performance of CARS. An 
analysis of the uncertainty in measuring the rotational tempera- 
ture is presented for diatomic gases. Inverse bremsstrahlung 
in plasmas is investigated as a possible process by which the 
parameters measured by CARS might be altered. An experiment 
is discussed in which an attempt was made to remove the back- 
ground signal inherent in most CARS experiments. Partial 
interferometric cancellation of two CARS signals was observed 
by the proper positioning of two sample cells in an otherwise 
standard CARS experiment. 
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I. 



INTRODUCTION 



Coherent Anti-Stokes Raman Spectroscopy (CARS) is a method 
of chemical analysis yielding molecular spectroscopic informa- 
tion by utilizing nonlinear scattering phenomena. As a part 
of the growing field of nonlinear optics, CARS has been used 
to measure nonlinear coefficients in solids, liquids, 
and gases [1-14]. Much of today's efforts are focused 
on the use of CARS as a nonreactive probe in environments not 
conducive to other techniques for spectral analysis [14] . 

The intent of this thesis is to assess the usefulness of 
CARS as an analytical tool using the example of diatomic gases . 
CARS is related to the more basic phenomenon of Raman Scatter- 
ing (RS) . Raman scattering involves a shift in the frequency 
of a small portion of the photons incident on a Raman active 
medium. RS has been utilized extensively to obtain molecular 
spectroscopic information. CARS yields the same type informa- 
tion as RS except that the collimated CARS signal is much 
easier to detect. 

CARS was feasible for chemical analysis only after the ad- 
vent of tunable lasers of moderately high power. The initial 
work with CARS was performed by Maker and Terhune [15] in 1963, 
but only recently has the field begun to expand rapidly, as 
evidenced by the large number of published articles. To per- 
form CARS in gases requires lasers of much higher power than in 
liquids and solids. Typically, pulsed lasers on the order of 
megawatts are used. 
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The use of such high power pulsed lasers has an inherent 
problem in that the power output from pulse to pulse is not 
constant. This introduces a source of noise that hinders the 
observation of weak resonances. In addition, in a mixture of 
gases, a CARS signal from other than the gas of interest is 
present as a background signal. These two factors represent 
at present a source of interference in observing weak molecular 
resonances . 

As an analytical tool CARS can yield information about 
molecular identity, concentration and rotational temperature. 
Numerical estimates for the performance capabilities of CARS 
are based on the ability to calculate the third-order nonlinear 
dielectric susceptibility, This susceptibility is com- 

posed of a resonant term, v_, and a nonresonant term, Xmd * 

Xo is calculable from the physical equations that define it 
but Xjjj^ was deduced from experimental observations. 

One source of observations for Xxtt> relied on results from 

^NR 

an experiment involving Stimulated Raman Scattering (SRS) . 
Appendix A contains a development of the equations related to 
SRS and shows the relationship of the SRS gain equation to 
CARS. 

A particularly difficult environment to perform convention- 
al spectral analysis with Raman scattering is in an ionized 
gas called a plasma. CARS is being used in plasmas and is 
producing better results than conventional Raman scattering 
techniques. An analysis is made to predict the extent of in- 
teraction between a plasma and the laser used for a CARS 
experiment in the plasma. In particular, the inverse 
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bremsstrahlung process for the transfer of energy is ana- 
lyzed. 

An experimental attempt was made to improve the basic 
CARS experiment by the elimination of the background signal. 
A description of the experiment and the results are included 
in this study. 
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II. 



NATURE OF THE PROBLEM 



The intent of this study was to theoretically investigate 
the analytical capabilities of CARS. The only area investi- 
gated was concerned with diatomic gases. It was desirable to 
estimate how small a concentration of gas could be detected 
and with what accuracy the rotational temperature could be 
measured. CARS has demonstrated its usefulness as a tool for 
spectral analysis. It offers orders of magnitude greater 
sensitivity than Raman scattering for low pressure gases. [16] 

Typically, CARS experiments are performed by focusing two 
pulsed laser beams of different frequencies into a sample cell. 
The beams cross at an angle that corresponds to the momentum 
conservation requirement. For gases this angle is sufficient- 
ly small that the experiment may be performed with the beams 
colinear. One of the lasers is tunable so that the difference 
in the frequency of the two may be changed. When the differ- 
ence in the frequency becomes close to an appropriate natural 
molecular frequency, the CARS signal increases markedly in 
amplitude . 

The intensity of the CARS signal is dependent on the pro- 
duct of the square of the intensity of one beam and the first 
power of the intensity of the other. Any amplitude instability 
in the beams leads to amplitude instability in the CARS signal. 
Also, the presence of an undesirable background signal masks 
very small signals of interest. 
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A brief sketch of several publications in the field of 
CARS has been included in this section of the study. The 
basic equations related to CARS have also been included to 
develop the framework on which the numerical calculations 
were based. 

A. BACKGROUND 

Maker and Terhune [15], in 1963, were the first to perform 

a CARS experiment. Prior to the development of good tunable 

laser sources most of the work with CARS was developed in order 

to measure nonlinear coefficients of materials. 

Yablonovitch, Bloembergen and Wynne [17] were among the 

first to use a tunable laser to generate a continuously tunable 

difference frequency. They investigated solid state materials 

such as InSb and GaAs . Levenson, Flytzanis and Bloembergen 

[1] followed closely with an experiment in diamond to measure 

(3) 

the variations in X as a function of the difference fre- 

quency. DeMartini [7] utilized CARS with H 2 gas in 1972 and 
Taran [4] in 1973 reported that he had used CARS as a means of 
determining the spatial distribution of H 2 gas in a flame. 

Taran 's recent work has resulted in the spectra of several 
gases, temperature measurements and spatial distributions of 
gases in flames [14] . 

B. NONLINEAR EFFECTS IN ISOTROPIC MEDIA 

Due to the nonlinear dielectric property of materials two 
laser beams may be mixed in a medium to generate a third co- 
herent beam such that 
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(1) 



0)2= 2 o)j^ - 0)2 

where o)^ is the anti-Stokes frequency generated from the 
pump and Stokes frequency and 0)2 respectively). 

The polarization vector associated with a material may be 
expressed as a power series: 



P (O)) 




E (o)) + 



X 



( 2 ) ^2 
(O)) (O)) 



+ X 



(3) -3 

(0)) ^(O)) 



+ . . . 



( 2 ) 



where X^^^ is the dielectric susceptibility tensor of rank 

i+1 and E is the applied electric field. The second order 
(2 ) 

term x is equal to zero for isotropic liquids and gases 

because of inversion symmetry. The lowest order term that is 

(3) 

present in a gas is therefore x / the third-order nonlinear 
susceptibility. There are several physical processes that in- 
volve the third order term, e.g. , 2o)^ + 

2o)j^ - ®tc. The last process, 2o)j^ “ ^2' pro- 

cess referred to as CARS. 

By assuming all fields are along one axis the components 
of the tensor, Xt and the electric field vector may be 
treated as scalars. The magnitude of the electric field may 
be expressed as : 



i (k . z 



Ei(o)) = 2 



- 0). t) 
1 



+ c. c. ] 



(3) 



where is the amplitude, k^ is the propagation constant 

(equal to m^nj^/c) , 2 is the length along the z axis and 

t is time in seconds. For the CARS process the term E^ in 



equation (2) includes 3E^ 
associated with CARS. 



E 2 [18]. The other terms are not 
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The resulting third order polarization, P 



(3) . 

IS given 



by ; 



p(3)_l 

“8 



(3) 

3x 



2 * 
01 02 e 



i [ ( 2 k^-k 2 ) z- ( 2 w^-o) 2 ) t] 



+C.C. 



( 4 ) 



The factor of 3 has been a point of inconsistency among au- 

(3) 

thors in the past. The notation for y (~w /W, ,o) ,a>,) is 
^ ^ a b c d 

consistent with that of Bloembergen [19] for the process 



0) 

a 



+ w + 



w 



(5) 



Solving Maxwell's equations for plane waves with polariza- 

(3) . . . 

tion P the anti-Stokes gain equation is : 




iiTto^ 

2ciTJ 





i {2k^ -k -k_) z 
)e 12 3 



( 6 ) 



2 

In terras of intensity, where | , after. integrating 

equation (6) over a length L, the result is 



^3 = 



y, 2 

4tt CO. 



V 



3X 



(3) 



2 2 2 
"2^ 



AkL/z 



(7) 



where Ak = 2k,-k.--k . The conversion efficiency is defined 

P 3 

as and if Ak = 0 the conversion efficiency becomes a 

2 

maxiimom. Because of the small dispersion of light in gas, 
CARS may be performed with , ^2 and colinear. For a 

focused beam (assumed to be a cylinder of plane waves and 
Ak = 0) the conversion efficiency becomes: 
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The third order nonlinear dielectric susceptibility, 
is composed of a frequency dependent resonant part (Xj^) ^ 

nearly frequency independent part, the nonresonant susceptibili- 
tyr 

X “ ^NR' 



X„ can be expressed [18] as: 
R. 



y = 2NC 

^ 4 

3noi- 



dn 



0) 



V 



[uj 



V 






- iy (a)^-oi2) 1 



( 10 ) 



When o)^-u )2 



0 )^, the resonant susceptibility becomes 



X 



R 



2Nc^ da i 
3*10.2“ 



(11) 



where N is the number density of molecules, 

Raman cross section for scattering, A is the difference in 
probability of finding a molecule in the ground and first 
vibrational states, and y is the line width, full width at 
half-maximum intensity. 
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III. NUMERICAL ESTIMATION OF 

(3) 

As defined xn Sectxon II. B., x is composed of two 

parts, the resonant and nonresonant susceptibility. The 
resonant part is associated with the various Raman active 
frequencies that characterize a molecule. The equation for 
X„ has been derived [18] from considerations of a driven 
harmonic oscillator and quantum mechanics. Xjgj^ is associated 
with the electronic structure of the molecule and is nearly 
frequency independent. 

Equation (10) defines Xc i^ terms of an entire vibra- 

X\ 

tional band, namely the ground state. To estimate y-d 

R. 

each rotational level requires an analysis of each of the 

\ 

terms in equation (10) that are functions of the rotational 
quantum number J. As will be shown in Section IV, the signal- 
to-noise ratio for a CARS experiment is a function of the peak 
value of Xj^* Being able to predict signal-to-noise ratios is 
important when deciding which analytical technique to employ 
for a certain situation. Due to the complexity of the calcu- 
lations a computer program was utilized to predict the peak 
value of X for several gases. 

X\ 

A. RESONANT SUSCEPTIBILITY, Xo 

R 

1. Cross Section 

The Raman cross section for scattering is typically 
quoted as the cross section for an entire vibrational band. 

For computational purposes it was necessary to apportion this 
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total value among the rotational levels . The relative popula- 
tion of a rotational level was used as the apportioning factor 
such that 

[ ^ (2J+1) exp[-BJ(J+l)/kT] / 

'^J Q \ dfi 

where 

Q = Z g_ (2J+1) exp [-BJ (J+l)/kT] , 

J 

J is the rotational quantum nximber, g is the nuclear spin 

u 

weighting factor and B is the rotational constant. The values 
used for ^ were from Chang [20] . His measurements of cross 
section were the results of experiments using a ruby laser at 
5140A. The cross section has a fourth power frequency depend- 

4 

ence (w ) thus an accurate calculation of X|^ requires that 
the data be corrected for this factor. 

2 . Linewidths 

On resonance the CARS signal is proportional to the 
inverse square of the linewidth. Two general extremes of 
linewidth were used, the Doppler broadened and collision 
broadened cases. The Doppler broadened lines are narrower 
than the collision (or pressure) broadened lines. 

3 . Rotational Transition Frequencies 

For molecules that exhibit Raman active vibrational 
and rotational frequencies the CARS signal increases in ampli- 
tude as the difference frequency approaches one of the molecu- 
lar resonance frequencies. The rotational frequency for a 
vibration- rotation transition of a diatomic molecule is ex- 
pressed as: 



) 



( 12 ) 



Total 
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( 13 ) 



0 ) T = 03 - 2o 3 X (v + 1) - a J (J + 1) 

v,J e e e e 

where 03^ is the harmonic vibrational angular frequency , 
is the enharmonic constant, v and J are the quantum numbers 
of vibration and rotation respectively, and is the rota- 

tion-vibration interaction constant. 

4 . Nuclear Spin Weighting Factor 

Homonuclear diatomic gas molecules have a statistical 
correction factor for the population of rotational levels due 
to the spin of each of the nuclei. Alternating lines in the 
rotational spectra have intensity ratios 1; (I/I+l) where I 
is the spin of each nucleus. The strong line of 0^/02 and 
N 2 are the even numbered J levels while H 2 ' s strong lines 
are associated with odd numbered J levels. TABLE 1 summar- 
izes the degeneracy weighting factor, g 

u • 



TABLE 1 



Gas 



NUCLEAR SPIN WEIGHTING FACTORS (g_) 

\J 



Nuclear 

Spin 



(I) 



9j 

Even # J Levels 



9j 



Odd # J Levels 



°2- 



1/2 



in_ 



^2- 



1/2 



Note that CO is not a homonuclear molecule and thus has no 
degeneracy factor; i.e. , g, = 1 all J. 

u 

5 . Computer Program 

To calculate Xd ^ nearly continuous function of 

K. 

the difference frequency required a computer program. A 
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Fortran program was written (See page 54) to calculate the 
maximum amplitude of Xj^ ^ function of the difference 
frequency according to the following equation: 



,2Nc da A 



R, (o)^ ^ 2 ) 



dQ Q 



y ; y , »T._ 

“v,J ■ -iY(“i-“ 2 ) 



(14) 



Tables 2 and 3 give the results of the computations. 



B. NONRESONANT SUSCEPTIBILITY 

'^NR 

In 1967 Rado [21] performed a series of experiments to 
measure the magnitude of the nonresonant susceptibility, 
of several gases. The impetus for the experiments was to mea- 
sure the peak resonant anti-Stokes^ power from a hydrogen cell 
and compare that to the nonresonant anti-Stokes power generated 
in a sample cell. The ratio of powers would be equal to the 
square of the ratio of the susceptibilities; 



sample ^ 



X (Sample) y 
NR J 

x'V 

P 



(15) 



The value Rado chose to use for Xp from an earlier 

work with phonon lifetimes and SRS gain coefficients. The 
relationship of Xp the SRS gain coefficient is discussed 



in Appendix A. His calculated value of 
.3 



'^NR 



for N 2 is 



1.35x10 . Based on Table 3, Xt^ N_ at 300 K 

erg , ^P 2 



-16 cm' 



is 2.4x10 

erg 

susceptibilities is 



Thus , the ratio of the peak to nonresonant 
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TABLE 2 



Gas 






^2 



«2 



CO 



Linewidth (cm 
. 0034 
.0076 

.0055 

.012 



. 036 



.07 



. 005 



.07 



Temp(°K) Xp(cm /erg) 



300 



8.3x10 



-15 



1500 



2. 0 X 10 



-16 



300 



3.3x10 



-15 



1500 



1 J. X 10 



-16 



300 



1. 0 X 10 



-14 



1500 



4. 3 X 10 



-16 



300 



3.3 X 10 



-15 



1500 



2. 2 X 10 



-17 



19 



- ♦ 




TABLE 3 



Gas 


Linewidth (cm 


Temp ( ® K ) 


Xp (cm^/erg) 




.16 


300 


2.6 


K 10-1« 


°2 


.07 


1500 


2.5 


X 10-1’ 




.14 


300 


2.4 


X 10-1« 


^2 


.06 


1500 


2.6 


X 10-1’ 




.03 


300 


1.2 


X lO"^"* 


«2 


.08 


1500 


4.0 


in-16 

X 10 




.15 


300 


2.6 


X 10-1« 


CO 


.011 


1500 


1.1 


X 10-1^ 



20 



178. 




An independent method to estimate the magnitude of 
was provided by the spectra of published by Taran [14] . 

For the unresolved spectra the difference in frequency of the 
maximum and minimum values of the CARS signal, 6o), is given 
by the relation 



6o) = y( 1 + 



Xp2 \2 



(16) 



4X 



Assuming a linewidth of 1.8 cm 
the calculated result is that 



NR 

-1 



and a 6o) of 42 cm 



-1 




47. 



(17) 



By estimating the level of the background signal in the 
N 2 spectrum, another ratio of equal to 




N. 



75. 



(18) 



Knowing the peak value of Xt, yields an estimate of 

R NR 

An earlier work by Taran [4] provides another means of 
estimating the ratio of the susceptibilities. The total 
pressure of a H 2 ,N 2 mixture was maintained at one atmosphere. 
The anti-Stokes intensity generated was plotted as a fiinction 
of the H 2 concentration. 
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scale linearly with pressure, then the 



If Xr and Xnr 
total susceptibility of 



the mixture may be approximated by 



x<^’ = XrP + [1-PIXnr 

where P is the partial pressure of H 2 . As a ratio of Xj^ 

to Y the observed anti-Stokes power is approximated by 
'^NR 



I 



3 






( 20 ) 



When the log of intensity is plotted against the log of P 

(See Figure 1) , the predicted and observed curves are in fair 

4 

agreement for a ratio of 10 :1. That is 



X 



(H^) 





( 21 ) 



Based on 
of Xp for 
estimated to 
ratio is 



the analysis of the H^/N^ plot and using the value 
H- from Table 3, the value of x,.tt^ for N_ is 

2 2 2 

be 1.2x10 Thus, another estimate of the 

erg 



Xp 2 4xl0~^^ 

P = — -- _ip 200. 

^NR 1.2x10 

The factor of four between the high and low estimates of the 
ratio is not surprising based on the variety of sources used 
to estimate them. 
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The solid 



Figure 1 



Taran's results follow dotted line 

X 

line is the predicted curve for — 

3 

two curves are equal from 10 to 



R 

NR 



= lo"^. The 
10^ ppm. 




H 2 Concentration (ppm) 



IV. STATISTICAL CONSIDERATIONS 



The anti-Stokes power resulting from a tightly focused 
CARS experiment is expressed as [18] : 



P3=K 



3X 



(3) 



2 T. 2„ 

P P 
1 2 



( 22 ) 



where 



K 



■fe) 



2 4tt 03- 






(23) 



The resonant susceptibility may be expressed in terms of its 

peak value, x / t»y; 

P 



= 

^R 2Ao 3 - iy 



(24) 



where Ao3 = 03, -03.,-03 . 

1 2 V 

For weak signals in which x absolute value 



squared of x 



(3) 



is given by 



X 



(3) 



X-nX 



= lx 



NR 2Aoj-iy 



(3) 



%R 



2xpYAo3 

2 2 ■ 
4 (Ao3) +y 



(25) 



For a 



The magnitudes of x'^^ its extrema are Xjgj^jlXp/2’ 

signal strength proportional to the difference in the peak and 

background susceptibility 



Xt. 2 



^S“^^NR '** 2 ^ ~ ^NR ^NR^P 



(26) 



Thus, for x-D << Xmtd' anti-Stokes power above that due 

IT ^NR. 

to Xj^R alone, is 



24 



(27) 



^ ^NR ^1 



where rirj, is an efficiency coefficient to take into account 
the optical losses, detector inefficiencies and observed con- 
version efficiency factor. 

Typically, the largest source of noise when performing 
CARS on samples at one atmosphere pressure is associated with 
the power fluctuations of the pump and Stokes beams. If the 
power of the laser source, , has a relative standard devia- 
tion of e , the standard deviation of the power is e P^ • The 

P P L 

fluctuation of the background signal power P^, due to 
is thus : 



'’n “ ^ ’‘nr 



(28) 



assuming each beam has the same relative standard deviation e^. 

For a low intensity CARS signal in the absence of back- 
ground one finds the major source of noise to be the shot noise 
of the signal. For an N photon signal the standard deviation 
is /W. 

For a pulsed signal the noise power due to shot noise is 



2 ^ 3 ^^ 



N 



(29) 



The signal-to-noise ratio for the background limited case 
given by equations (27) and (28) becomes 



S/N 






(30) 



whereas the signal-to-noise ratio for a situation in which 
there is only a very weak signal is given by equations (27) 
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and (29) to be 



1 1 

S/N (K Xp^ (T/hv)^. (31) 

As an application of the signal-to-noise problem consider 
the measurement of Xp from a spectral envelope. The uncer- 
tainty of Xp, Op, is expressed as; [22] 



a ^ ^ a. ^ 

P 1 



(32) 



where is the standard deviation in determining tp from 

th 

the i pulse of a pulsed laser. For CARS with a background 

limited situation and uncertainty in the signal 

amplitude is approximately a constant, associated with 

the laser power fluctuations. 

3Xp 

Thus ( ) where a is basically the 



m 






m 



(w. J,T) 
1 / 



noise power given by equation (28): 



a = 3 e K x^jT, 
m p '^NR 1 2 



(33) 



I, T is the intensity of the spectral envelope at co. , 

( (l) , J f 1 ) 1 

1 / 

due to the quantum rotational number J, at temperature T. 
The result of this analysis yields [16] : 



a 



a = 



m 






/ 2Agj > 2 
' TTY 



(34) 



or 

1 

_ 1 .2Ao). 2 

^p S/N ^ TTY 



(35) 
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In general any standard deviation is dependent on the nuniber 
of times, N, the measurement is taken such that 



a 

P 



1 



a — 






The conclusion is that for a given S/N the term is the 

effective number of pulses required at intervals of Ao) to pro 
duce an uncertainty, a^, in the measurement of Xp- 

The uncertainty in determining the rotational temperature 
is a product of the nearly constant background uncertainty, 
and the small change of intensity with temperature. A 
convenient definition for a is: 



m 



a 



m 



I (5, 
max I 



(36) 



max 



where 6^ is the relative uncertainty in the amplitude of 
■^max 

the line of maximiim intensity. The uncertainty in the tempera 
ture then becomes [16] : 



1 






(37) 



The intensity of a line is proportional to the Boltzmann 
population distribution. The change of intensity with tempera 
ture has been shown [16] to be; 



^^J,T ^ [BJ(J + 1) - kT] 

3T J,T ^^2 

The resultant standard deviation in temperature is 
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It has been estimated [16] that for an uncertainty in the 
maximum intensity equal to 2% the temperature of CO at 1500®K 
may be determined to within 10 °K. 

This uncertainty in temperature is for an unrealistic 
situation in that each rotational line is completely re- 
solved. A' more realistic situation is the case of 
a continuous envelope as given in equation (14) . A computer 
program (see page 59) was written to calculate the uncertainty 
in a temperature measurement, based on equation (37). 



, 3I(o). J,T) 2 

y ( i_£ ) 

^ 2 i ^ 3T ^ 

m 



o = Ie = Ky_£ 
m p p ''•P p 



, 3(K J,T) 




For the same conditions used to calculate the previous 
uncertainty in temperature the results of the envelope calcu- 
lation yields an uncertainty of 15 °K. This value is of the 
same order of magnitude as the statistical prediction of 10°K, 
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estimated by using equation (39) . It is to be expected that 
the uncertainty would be higher for measurements on the 
spectral envelope. Verification of the statistical approach 
means that the much simpler method for predicting the uncer- 
tainty is a good approximation. 
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■ V. CARS IN' A PLASMA 



Low pressure ionized gases, or plasmas, represent a class 
of materials that are receiving very close investigation. The 
population of the various energy levels reveals much about the 
plasma. Determining the population of levels is a very diffi- 
cult task with conventional Raman spectroscopy. CARS has been 
utilized to obtain spectra of several plasmas at the Naval Re- 
search Laboratory, Washington, D. C. The experience there 
indicates that the technique is not difficult and produces 
useful information. 

Of importance when analyzing CARS in a plasma is how much 
the high power laser beams may_ be expected to perturb the 
system under investigation. The major source of interaction 
between the plasma and laser beams appears to be from the in- 
verse bremsstrahlung process [23] . The free electrons in the 
plasma will gain kinetic energy by this process. 

If the laser frequency is above the plasma frequency the 
light will propagate through but there will be some attenuation 
or absorption. The plasma frequency is given by [24] 



1 




-3 

where ri^ is the number density of electrons (cm ) , e 
is the charge of an electron (esu) and is the mass of an 

electron (grams) . 

An estimate of the change in the electron kinetic energy 
may be made by assuming there is an absorption coefficient, a. 
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for the plasma. The energy absorbed per unit volume of the 
plasma is 



"ABS 



T- ,, -ox, ^ 
I (1 - e ) At A 
o 

Ax 



P^Ata 

o 



(41) 



where I is the incident intensity, P is the incident 
o . o 

power. At is the pulse duration, A is the beam cross sec- 
tion and X is the path length of the medium. If all the 
absorbed energy is transferred to the kinetic energy of the 
electrons then the energy absorbed per electron is 



aP At 

E — ^ 

ABS An 

e 

As a fraction of the initial energy, 3/2kT^, 



^ABS 

3/2kT^ 



aP At 
o 



3/2An kT 

' e e 



(42) 



(43) 



where is the electron temperature. The 'absorption co- 

efficient, a, for inverse bremsstrahlung has been reported 
as [25] 

1 ^2 6 

® ^ ^3/2kT^ r 3/2 3 “ exp (-j^) ] (44) 

h c M ' V 
e 

where n^ is the ion number density, z is the charge of the 
ions and v is the frequency of the incident light. Evalua- 
tion of equation (43) with a defined by equation (44) results 
in a relative change of temperature equal to 

3T^ 1.3xl0^^n^z^P^At[l - exp(-hv/kT)] 



T 



3 



3/2 



( 45 ) 
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This agrees favorably with an equation for the change of 
electron temperature given by [23] 



AkT 

kT 



25 2 

1.4x10 n-z [1 - exp (-hv/kT) ]/ P^dt 
e 



( 46 ) 



where electron volts and is the power of 

the laser in watts. 

The equipartion time, t^, is the average time required 
for the electrons to transfer their excess kinetic energy to 
the ions and thus reach an equilibrium temperature. The equa- 



tion for t^ is given by [25] 



t = 
e 



25.2 A(T ) 
e 






3/2 



(47) 



where A is the atomic weight of the ions , T^ is the electron 
temperature in degrees Kelvin, and is in cm If the 

duration of the laser pulse is shorter than t^ then any 
energy the electrons gained by inverse bremsstrahlung will 
not be transferred to the ions and there will be very little 
perturbation of the molecules under investigation. 

Consider the following as a typical example: 

4 

plasma; Deuterium at 50 Torr,Te =10 K 

11 -3 

7e = 10 cm , 20 cm long 
• -4 2 

laser; 5500A, 5MW, 10 cm beam area, 20 nsec 



duration 

The plasma frequency, 

The equipartion time , t^ , 



equals 1.8 x 10^^ sec ^ and • 

L p 

equals 20 micro seconds. The 
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relative change in electron temperature given by equation (46) 
is .0001%. Little perturbation of the plasma is expected from 
these conditions. 
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VI. BACKGROUND SUPPRESSION EXPERIMENT 



As discussed earlier, the nonresonant susceptibility, Xmd' 
is responsible for the background signal. The power amplitude 
fluctuations of the lasers produce noise on this signal. The 
resonant signal above the background also has noise on it. It 
would be much easier to detect the resonant signal if there 
were no noise on the background or, better yet, no background 
at all. The desire for reduction of the background signal led 
to an idea [26] for an experiment to determine if it was possi- 
ble to completely suppress the background. In addition to 
yielding much better signal quality the elimination of the 
background would limit the detectibility of CARS to the limit 
of conversion efficiency. 

Typically, the diluent of a sample provides the largest 
contribution to The basic idea for suppression of the 

background signal was to perform two CARS experiments in series . 
In the first experiment the generated anti-Stokes beam, co^/ 
would be due to only ^ diluent) and in the other 

I 

experiment the anti-Stokes beam, would be due to 

Xj^ (the sample plus the diluent) . If the to^ beam, from the 

first cell, were then forced out of phase by Ntt radians 

I 

(where N is an odd integer) relative to the second beam, 
and superimposed on the second then the resultant signal should 
be only that due to Xo sample. 

Molectron Corporation of Sunnyvale, California, through 
Dr. G. Klauminzer provided equipment, facilities and assistance 

for conducting the experiments . 
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Figure 2. Arrangement of Equipment for Background Suppression 
Experiment. 

Figure 2 shows the basic arrangement found to be most 
convenient. A pulsed nitrogen laser was used to drive two 
tunable dye lasers in an arrangement reported by Dr. 

Klauminzer in Reference 12. 

Benzene was chosen as the first material to work with 
since it was known that the CARS signal was very strong and 
visible. For benzene the values for W 2 and were 

480., 504. and 458.2 nm, respectively. 

Both cells were filled with benzene in anticipation of 

f 

demonstrating complete cancellation of with co^- It was 
known that a translation of cell one along the focal axis 
would change the intensity of by changing the volume in 

the focal region but it was also discovered that the same 

I 

translation changed the phase of relative to To 

f 

completely cancel it was known that several conditions 

had to be met. The beam had to be ir radians out of phase, the 
intensities had to be equal and the beams had to be coinci- 
dent. A plot of the resultant intensity versus the translation 
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of cell one revealed that there was constructive and destruc- 
tive interference of the two beams. The destructive interfer- 
ence was not complete though. A visual observation of the 

I 

co^ + beam revealed interference fringes . By steering the 

I 

beam off of the beam the fringe pattern could be 

changed. It was suspected that chromatic aberration in the 
focal lens was responsible. Calculation of the effect of 
chromatic aberration indicated that the fringe pattern should 
be approximately that observed. The use of an achromat for 
the focal lens reduced the number of fringes visible but it 
was found that a slight rotation of the lens changed the 
pattern considerably. 

The best measured ratio of minimum to maximum intensity 
as a function of cell translation was observed to be 10:1. 
Since it was the electric fields that were interfering, the 
10:1 ratio amounts to about a 60% cancellation factor. The 
cancellation factor was arrived at by considering the incom- 
plete phase cancellation of the electric fields. The observed 
intensity is proportional to 



i0. 



i0. 



E. 



+ E, 



(48) 



If the amplitudes E^^ and are equal, the resulting 

intensity is proportional to 

P ie i0„ -i0, -i0„ 

E (e -^ + e (e -^ + e (49) 



or 



a 



2 

E (2 



+ 2 cos (0, 



- 02 )) 



(50) 
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results 



Complete constructive interference, when 0^^ = Q^, 
in an intensity four times that of one beam, and complete de- 
structive interference would result in no intensity. 

It is not surprising that on the average complete inter- 
ference does not take place. For a small difference in 0^ 

2 

and 02 the observed magnitude is slightly less than 4E or 
slightly more than zero. 

An observed 9:1 intensity ratio is equivalent to a 3.6 to 

2 

.4 ratio on the normalized scale, E =1. Thus, the inter- 
ference produces a signal equal to 40% of what one beam would 
produce, therefore, a cancellation factor of 60%. 

Next, a three molar sodium benzoate solution was used in 
the cells in place of benzene. Again there was incomplete 
cancellation. The best observed ratio of intensities was 5:1. 

It had been predicted that sweeping the frequency of 0^2 
record a spectriom would not alter the phase difference of 

I 

and The first observations made while sweeping 0)2 re- 

I 

vealed that and went in and out of phase over a 

period of about 10 cm By changing the optical path length 

of 0)^ or and 0)2 by use of microscope slides the 

period could be increased to about 60 cm 

The observed increase in the period of the oscillation 
verified that the optical path length could be controlled to 

I 

shift 0)^ and in and out of phase. With the appropri- 

ate optical path length the period of the oscillation could 
be adjusted to infinity. 

A spectrum was recorded for the three molar sodium benzoate, 
two cell, series cancellation experiment. Of interest from the 
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spectra is the interference of the background with respect to 
the noise levels observed. Figure 3 is a portion of that 
spectrum depicting the observed resultant signal. When the 
interference was constructive the noise was much higher than 
for the destructive interference. 

As stated earlier, the reason for the experiment was to 
increase the sensitivity of CARS. An analysis of the signal- 
to-noise ratio for a cancellation experiment indicates that 
the cancellation has to be much more effective than was ob- 
served in this experiment to gain an order of magnitude in- 
crease in the signal- to-noise ratio. 

2 

Assxime b is the fraction of background intensity remain- 
ing after cancellation. Then, the S/N ratio is given by 



S/N = 






t 



similar to the development of equation (30) . An increase in 

2 

S/N of 10 requires b = .01, or a 99% cancellation factor. 

A 400:1 intensity ratio would correspond to such a cancella- 
tion factor. 

The experiment as performed was very sensitive to the 
alignment and quality of the optics involved. One of the 
major contributions to the incomplete cancellation is probably 
the diffraction limited ability of the center lens to perfect- 

I 

ly image the initial tu^ beam onto Other macroscopic 

factors such as turbulence and inhomogeneity of the samples 
are other obstacles to complete cancellation. 
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background signal 



VII. RESULTS AND DISCUSSION 



The basis for determining the areas of application for a 
new analytical tool, such as CARS, must depend on recognizing 
the limits beyond which the technique cannot be expected to 
yield the desired information. 

The results of this analysis are presented as lower limits 
of concentration that are expected to be detectable by CARS. 
For a mixture of gases, the detectability of one gas depends 
on the ratio of Xo Xmt,/ the same as Taran's plot for 

H 2 and N 2 « Plots of the ratios for other gases with N 2 as 
a diluent are presented in Figure 4. The minimxim detectable 
concentration is dependent on the required signal-to-noise 
ratio for a given experimental condition. The lower limit of 
detection for a pure gas at reduced pressures would be limited 
to the conversion efficiency if not for the photon statistics 
involved. From the signal-to-noise analysis, estimates of the 
detectable lower pressures are summarized in Figures 5 and 6. 

Based on the results of Section IV the rotational tempera- 
ture of a gas is expected to be measurable to within 1%. That 
assumes the amplitude fluctuations of the CARS signal observed 
remain less than 10%. 

In a plasma, the inverse bremsstrahlung process is ex- 
pected to be a relatively insignificant perturbation to the 
parameters determined by a CARS experiment. 

The experimental attempt at cancellation of the background 
confirmed the idea that interferometric cancellation is a 
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possible technique for improving CARS when the background be- 
comes a problem. The observed results were not of the magni- 
tude hoped for but did serve to verify the concept and to 
recognize a possible source of limiting performance. 

Turbulence in the samples and inhomogeneity may be part 
of the reason that incomplete cancellation occurs. The major 
contribution to the lack of complete cancellation is probably 
the result of diffraction limited optics. The first beam 
cannot be refocused through a lens and have precisely the same 
geometry that it had originally. 
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Relative Anti-Stokes Intensity 



I 



Figure 4 




Concentration [ppm] of Sample 



Log-lo^ plots of N_ gas mixtures at 1 
total pressure. at 300 K. Upper CO 

300 K, Lower coplot at 1500 K. of 

to be 1.2 X lO"^® at 300 K, and 

-IS 3 

.24 X 10 at 1500 K. 

erg 



atmosphere 
plot at 
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Figure 5. 



Signal-to-noise ratios for gas mixtures at 300 K 
with a total pressure of 1 atmosphere. Background 
pressure and signal due to N 2 * The area boundar- 
ies are defined by assuming values of the un- 



certainty in the pulsed laser intensity, 

-2 -1 
10 to 10 . 



£ , from 
P 
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Signal-to-Noise Ratio 




PRESSURE (ATMOSPHERES) 



Figure 6. Signal~to~noise ratios for Doppler broadened 

gas at 300®K as a function of pressure. Noise 
due to photon statistics. The assumed parameters 
are 

= .01, P^ = 5MtV, P2 = .5MW, T = 20 nsec, 

X = 5500 A. 
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VIII. CONCLUSIONS 



This analysis presents a model by which the limits of the 
information available from CARS may be predicted. At STP 
typical diatomic gases should have a signal-to-noise ratio of 

-3 

1 at partial pressures of about 10 atmosphere. Hydrogen, be 

cause of its narrow linewidth and population distribution, 

“5 

should have a signal-to-noise ratio of 1 at about 10 atmos- 
phere. Pure gases at low pressure where the linewidth is 
determined primarily by Doppler broadening, may be detected 
at pressures several orders of magnitude lower than this. 

The ultimate detectability depends upon laser power and satura 
tion effects. 

Simple relationships have been derived which allow the 
prediction of the uncertainty in temperature of a gas for an 
assumed signal-to-noise ratio. 

In a plasma the inverse bremsstrahlung process is found 
to transfer insufficient energy to the plasma to significant- 
ly modify the parameters measured. 

The experimental attempt at background reduction by inter- 
ferometric cancellation verified that the principles were 
correct. Complete cancellation, however, was not observed. 
This is thought to be due to the diffraction limitations of 
optical lenses. 

One area of additional research and experimentation that 
could further extend the applications of CARS is that of the 
background suppression. Another area of improvement needs to 
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be that of the fluctuating intensities of the high power 
pulsed lasers that are available. 
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APPENDIX A 



STIMULATED RAMAN SCATTERING 



A. DEVELOPMENT 

The interaction of two electromagnetic fields in a material 
is such that different combinations of the two may produce 
additional fields. The two additional fields of interest here 
are the Stokes and anti-Stokes frequencies. 

Raman scattering has been observed extensively and found 
to produce the Stokes and anti-Stokes frequencies at plus and 
minus the vibrational frequency, of the molecule from 

the incident light, 03 ^. 



"'stokes - "'l ■ "V 



"'anti-Stokes "'l 



+ CO. 



V 



(Al) 

(A2) 



Stimulated Raman scattering predicts a gain in the strength 
of a field as it traverses a media. The media of interest in 
this study are gases. The gain will be related to the non- 
linear susceptibility of those gases. 

Beginning with Maxwell ' s equations , which contain various 
assumptions implied by their form, the plane wave equation is 
derived from the curl of the curl of the electric field, E. 
Assuming a solution to the wave equation of the form: 



, -i (kz 

E = I (e e 



cot) 

+ c. c. ) 



(A3) 



where e is the amplitude , k is the propagation constant 
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and w is the frequency. For simplicity, all fields have 
been reduced to one dimension along the z axis. From the 
nonlinear polarization term the various possible frequencies 
are defined and their field properties (gain or attenuation) 
are explained by computing the solution to the wave equation. 
For completeness both MKSA and CGS results are presented. 

The development in MKSA units begins with Maxwell's basic 
and constituent equations. 



V X H = J + D 



D = 



e E + P 
o 



V X E = -B 

V • E = 0 



J = aE 



P = X E + P 



NL 



B = 0 



B = 






The curl of H becomes 



V X H = aE + e^E + P = aE + eE + 

Then the wave equation, defined from VX(VxE) is 



-V 



dt 



(V • B) 



d 

"’^o dt 



(V X H) 



or 



E + 




V^E 

^o" 



£ NL* 



(A4) 



Using the assumed form of the solution, equation (A3) , and 

2 2 

the assumptions that a = 0 and 3 e/3z is ve3ry small, 
then equation (A4) becomes 
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i (kz - wt) 



0) , 
(ee 



2 2 

. „ ^ V , k c , i(kz - o)t) ^ . 

+ c.c.) + y- (ee ' + c.c.) 

2n 



kc' 



o Z 



. c.^ = 



- i i 

e NL' 



(A5) 



The usual definitions: 



2 1 
c = 



u e 
o o 



r) = /e/e 

o 

k = 0 ) ri/c 

have also been used. 

For the two interactive fields and the nonlinear 

polarization, P is expressed as: [18] 

NLj 



P = 3e 
NL o'^ ^ 



1 2 * 



(A6) 



Three frequencies are possible from the E^ E 2 combination; 
0 ) 2 , 2o3^-aj2 / 2tu^+a)2 • If E 2 is the Stokes field with frequency 
0)2 (assuming oj^< ^ 2 ) then a resultant field, E^, of fre- 
quency 2aj^-032 is the anti-Stokes field. For the anti-Stokes 
calculations P^^j^ is given by: 

o * -i [ (2k, -k_) z- (2 o),-u 3«) t] 

PnL ^ ^^o^ ^¥^1^1^2 c.c.)]. (A7) 



The 



resultant field may be E 2 for the case when 



3 ,1 



'’nL ' ^2' 



-i (k2Z - 032t) 



+ c.c. ] 



(A8) 
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1 . Stokes Gain Equation 



To solve equation (A5) for the Stokes condition, equa- 
tion (AS) , there are two possible alternatives that lead to a 
gain equation. One way is to assume that 9 £ 2/22 is very 
small and k 2 is complex such that 

k2 = k' + ik". 

The imaginary part of k 2 , k” , is then evaluated as: 

U) 3x"lE,°l^ 

V « — _£ i 

^ 2cn 

where x" is the imaginary part of is the 

magnitude of at z = 0. 

The observed intensity of E 2 is proportional to the 
2 

time average of IE 2 I or 
where 



K T. * 1/ ik'2 -k"z ^ * -k'z -k”z, , ^ , 

® ^2 ® ® (C.C.) 



Thus, 1^2 1 proportional to e . x” is a negative 

— 2k*'z 

quantity [27] and therefore e is an exponentially 

growing function with distance. 

The other possible way to solve equation (A5) is to 



let 



and then solve for 




9£2/9z* The result is that 
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2 

3(X' + ix") |Ei°| 



1 

2^2 



] 



and again for x" negative there is gain for increasing z. 

Both approaches result in an intensity gain coeffici- 
ent, g, such that 



I 



2 



I 



o gz 
2 ® 



where is the intensity of at z 

that 



g = 




|3X”I 



2 



0 . The result is 



(A9) 



For cgs units the gain is 



g 



477032 



H2C 




(AlO) 



2 . Anti-Stokes Gain Equation 

Solving equation (A5) for the anti-Stokes condition 
(A7) results in the gain equation for the anti-Stokes field. 
The only way to solve equation (A5) is to assume 




The result is that 



3^3 _ ^“3 

3z n 3 C 



3(X’ + ix”) 



1 * iAkz 

8 ^ 1 ^ 1^2 ® 



(All) 



where Ak = 2k^ - k 2 - k^. 

Again if x" is negative, there is gain. In cgs units 
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(A12) 



9z (cgs) 



3e . 



= 4tt (- 



9z 'MKSA. 



B. CORRELATION OF VARIOUS ARTICLES 

The only published values for nonresonant susceptibilities 
of gases are those given by Rado [21]. It was desirable to 
check for consistency in the equations used by Rado and the 
ones used to calculate resonant susceptibilities. 

The basis for Rado's experimental values of '’^he 

resonant value of ^2 other experimental 

data [23] . The gain coefficient for Stokes Stimulated Raman 
Scattering provided Rado with the needed value of Xn 
at pressures greater than 10 atmospheres. The gain coeffici- 
ent was given in [28] in cgs units as 



k [S tt^L a Na^] T 
g _ r 

I, e cn "K 



X 10 (cm/W) 



(A13) 



where L is the local field correction factor, N is the 
number density of atoms, A is the population ratio, is 

the phonon lifetime, a is the polarization coefficient and 
Ij^ is the incident laser intensity. 

An independent article by Maier [29] used the following 
definition for the gain coefficient 



g _ 8itN (Jpg 
^1 'hn^ Av 

where Av is a linewidth. Rado's gain is equivalent to Maier 's 
if 

A = 1. 
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L = 1. 



and 



e = 



n 



2 



T = l/2irc Av. 
r 

The experimentally observed gain for H 2 above 10 atmos- 
pheres was stated as 1.45 x 10 cm/W for* an incident light 
source of 975 nm. If the gain coefficient developed in equa- 
tion (AlO) is divided through by intensity then 





2 2 
n 



|3x"l 



Solving for 1 3x" | the result is 



(A15) 



gcA 

|3x"l = 4 (A16) 

2 

assuming n =1. The value for 3x” is the same as that 
used by Rado to represent x”* 

The gain coefficients given in two textbooks more nearly 
resemble the gain coefficient in equation (A9) . Pantell and 
Putoff [30] define g in MKS units as 



g 



47t x" 



(A17) 



whereas equation (9) may be rearranged as 



471 |3x"| 






(A18) 



differing from Pantell by a factor of If - 1/ they 
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become the same. Yariv [27] gives a voltage gain coefficient 



as 



g' 


= % |E 


1^ lx" 


1 9 


(A19) 


The intensity gain, g. 


would be 


2g'. 


Thus , the 


intensity 


gain coefficient is 










2g' = g = -^ 


|E|^ |X”I 


= 

n,c 


|X"I |E|^. 


(A20) 
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PROGRAM TO CALCULATE THE LOCAL PEAK VALUES OF TPE 
RESC^ANT THIRC-CRDER NONLINEAR CIELECTRIC 
SUSCEPTIBILITY, CHt FOR A DIATOMIC GAS 



















X — 
















X 


(X 2. 
















OJ 


X— o 


































fH 


ox*^ 


• 














W 


Z'^o^ 


UJ 














> 




o 














o 


















z 


l->-wx 
















UJ 


ox i 


tj 














zx 


<IZXX 






1 




U1 






CJO^ 


xxox 


h- O 




1 




H- « 






H-4X 


xoxo 






1 






' X 




»-QC — 


»-cya:z 


< 




1 




z 


— < 




ou-x 


ZXX>-4 






1 




o 


X o: 




X OJ 


►-IQC 


00 o 




1 




•-4 


OJ o 




00 


xoo>* 


k: CD 




1 




t— 


o 




<1.^ 


Z XO-J 


OJCXQC 




1 




<*. 


Z rH oc. 




i/JZ^ 


OJXZZ 




• 


CU 1 




-J 


l-« w 




orioi— 


l-f^*-4X 


ULLJ 


— • 


1 




-3 


-> 




UI-IZ 


>— ZH»0 


>* 


o 


UJ 1 






•J X X 


X 




<io o 


cr»-4Lj 




O 1 




*J 


X OJ X 


o 


u<t»-^Ha:>-4H-x 


<txz 




O 1 




< 


X »- 




QCOOXXJ— ZQl 




oc 


O 1 


•^cvirn^ 


o 


oox 


i/) 


ZXZX>-4<1<1X 


00 


UJ 


1 






00X00 oo 


»- 




UJ z 


H-4 


1 




< 


XQ^CC X 


oo 


X>00 1 OQOOX 


OJUULiJ 


X 


1 






XXX— a 


l-H 


<l 


ZZI-4ZO 


UJI-O 


H-4 


1 




oc 


a:x>oi- 


00 


(XCtX 


a>oz 


2:<^o 


•- 


1 




O 


0DXZ^V7 


z 


<t<Xi-i ox 


^o: 


z 






X 


xoo>-» 


o 


XXZOH-O oc 


UJ30J 


UJ 




UJ 




oo X • 


OJ 


<t:Do<t<L»-»ox 


xo> 


o 




o 


>- 


3izao 




|-0»-aX»-Z*-^X 


K-JX 


1-4 




H-4 




Z»— *-40 


<t 


OUJh- 


oaz>-» 








X 


X 


»-H<l ox 


H- 


H»x<ixcx.zoa 


000-rz 


z 


1 


z o 


z 


— o 


<0 


oi-a OCX 


<X UJ 


<t 


1 


UJ Z 


1-4 


— z<xcr 


o 


xzo 


x<txx 


XOO 




1 


zouzo 




I.H-I2.XX 




X QCZXX<X 


h-O 


00 


1 


UJOUJX 


Q 


w <h^q: 


X 


l—X ( 


H4»- XJ- 


X Qt: 


< 


t/) 1 OQCO 


UJ 


X— 0X0 


<l 


XxCL 44Z 




X 


< 1 




OC 




z 


• •^X(0 ••x^ •• 






o 1 


XI-4QC.LJ 


•-4 


0^X1- <l 


Qd 


o K 


XCJ o 


CJ—liC 




1 


cj<^ao3 




>-X|— zx 


X 


»-4 •• 


••XX ••z 


o<t ^ 


<c 


1 


>-a: 


OF 


<QlOXX 


»- 


oo> ••i-4xa:x»-4 


cr^z 


U) 


1 


X<I 


X 


cxo>-»xx 


z 


cjzxQt:<xx:s 


iX.(XLU 




1 




(X 


XOOJC*^^ 








ujo:j 


X 


1 






XCOX X 








UJh- >• 




1 




X 


ZXZOOI- 


X 






XZX 


< 


1 




X 


X(Xh-i<1 


X 








UJ 


1 




1- 


h-X— 


h- 







I 






55 



OOOOOOOOO 000«3000 
•^cvjro^u>'on-oou' (\jp>>TtnNur-co 

OOOOOOOOO 

ooooooooo ooooooo 

OOOOOOOOO OC7OOOO0 




OC3O0OO0O 

cn rn ’4- so u 
*^rg t\j c\j c\j CM (\j c\j 

oooooooo 

oooooooo 



0U?0<-J0U00 



oo 

P-CO 

(MCM 

oo 

oo 



oo 



oooooooo 

cj'0'^cMrn>r‘^^u 

oamnu'u'ifnmrt^ 

oooooooo 

oooooooo 



oooooooo 



o 

CM 



o 

o 

► • 





in ^ 

•. OCP 
*orci 

sO^ •CJ' 






P>0 »KM 




r-4 


•otn • 




xO 


O 




w 


00 •»04i-4 




-5 


moo • •» 




o 


•^ouoro 




•» 


•m-CCMCM 


— 




U' •p^OO 


V 


ni 


r->T •» • 


» • UJ 






— Qc. 


c/> 


l-H Pi^^sO 


• o 


a: 


or- ••-1 


XCOC/)— CO 


< 


•oorsj 


LUUJ« CO 


o 


•» #o ^ 


</)LUUJXULI 


•» 


•oo%rm 


aicro: oc 






ujcj>0(yiQ* 


0) 




OLOIUJUJ 




r*cNc»r* 


OOOUJ 


►— t 




— Ol otco 


jS 


in *^<10 


^TO ox 




>ror^ • 


t - • LUO 




<PT^O •• 


O •• ••Oh- 


O 


•f-4 • • 


oo < 


>o 


•Monicr 


CL • • 




•» • ^CT' 




-5 


cn •* 


• •u.x *o 


or 


o • *r- 


Qin •• «o • 


— UJ 


UJ«-» **04 


LULL* • «m 


OOL 


O ^>Tp-4 


CO •• mu- 


>ou- 


inr-ocM 


— O- h-h-x * 


^ •> 


inaDODCj> 


fO « ~ 


X— 


•o^ • • 


• XH- 


<j;o 


— <u>o^ 


OI-<2Z CX 


atvO 


<► tvO •• 


P^Q LULU o 




'4*^ •*<T| 


U-I-I200I-X 


••n 


m -.cno*^ •*3UJOO< 



—INI of^OOKMoouJOaro^ 

0 •• O •ujrn«-4 — 

:CsO— •o\n>T «Mjn»-iX»-i>-o» 
c»wo Mnmi-^o^-^-JOZXo 

kC>C ,-4 CMrOi-l • toi- - - - • 
tu«q; •o *04 •cMir\ •*•»•»•••••.•*•• 
asicT' •-» *^xx>c><i;<xxx 

1 i w (XNstn ii' !■ I 

X X a. 

aioOOOOOOf-»— h-F-h- ^-l— 

O • •<Kl<<<X<t<<l<l 

v^'tUMUC.LX.U.ULLt.CZ.ULULLX. 

• • «tJLJOOOOOOO 
OQCac Lj*-toi«^ii.i^ii.u.a.u.a.ULU. 
iNOisr 

ooooooooo 

ooo*^'^i'i>roo 

^rvjinu^u>um\Nur^ 



o 

X 

H-t 

cc 



oooo 







<l 




— 








X 




X 












o 








cx 




<w 








o 




X 












l-H 








o 




0^ 








< 




1— 








X 




< 












X 








o 












2 




<l 








H>4 




1- 








H- 




<l 








X 




o 








o 












1-^ 




UJ 








LU 




X 








3 




1— 








2 




X 






1— 


H>4 




o 






-iiP i*P 


X 




a. 




CO 


^ oooo 


</) 




X 




LU 


P» »-iCMO>Nr 










o 


<1 incnmin 


a: 




< 




ou 


X 


<L 




h— 




< 


<I O'^'O '-O 


UJ 




< 






O ww«mpwQ. 


o 




O 




a: 


* UJLUUJIUO 


o 








< 


Q. 1 ■! 1 ■ 


3 




Q 




> 


p> i-M n II icO 


2 




LU 


— 




H— ^cxQiQca: 




P-4 


a: 


— p — **o — 




— 0-SJC3CJ;— 


UJ 


nU 


HH 


O O •• O 


• o 


oo-^-^-^— • 


X 


•» 


3 


*r-*— •* 


Q. 


or-p^cMm-^o 


H* 




or 


i-iO— O 


2 


p^ ^ • t • • 




tl • 


X 


w p»<0 ••^X 


1— » 


••'Ooroooo 


UJ 


Of-i 


cr 


xm •inxooDx 




inwUJUJLULULU 


INI 


II 




o^>r^o It ^o 


LU 


WUJ • • • • • 


HH 


sO^ 


UJ 


11 X^X II XX II 


X 


LJh-ooooh- 


o 


O 


X 


00X0X000 


h- 


<CC l-H-* w ^^^O 


<!► 


O^ 


h- 


1-1 II O It XX II ^ 




«ajULU.CUU-U.I 4- It 


(p-« 


U-IO 




II »-ixxx»-« 


o 


0L,Jp> *1 M *1 I* 12 


►- 


o 


O 


O^XX<lXX^ 






HH 




o 




UJ 


o 


2 


>u 


-J 




LC 


o 


— 




X 




•ft- 


m 


* 




it* 




it* 




it* 




it* 




ooo 


OOO ' 


OOO 



56 



o 


0000^0 


oo 


ooooo 


o 


o 


ooo 


o 




GQo>o^r\irn 


^tn 


sOP-cno^ 


C\J 


O) 


'Tin^u 


r- 


cn 


mm>T'T>T'T 


>r>r 




in 


m 


unnm 


in 


o 


oooooo 


oo 


ooooo 


o 


o 


ooo 


o 


o 


oooooo 


oo 


ooooo 


o 


o 


ooo 


o 


l-H 


M il <1 1 


i-hhh 


i-Hi !>■<♦ il < 






1 <1 Hi 


•-H 


X 


xxxxxx 


XX 


XXXXX 


X 


X 


XXX 


X 


o 


oooooo 


OC) 


ououo 


o 


u 


oou 


<-> 



;> 

UJ 



ZL 

Ul 



























—1 




X 














h- 




















<l 




;> 






3 








O 










CL 










Z 




X 




H- 










oc 










X 










o 




X 




v» 




















O 










H -4 








CD 


X 




o 




X 










X 










1 — 




X 




•Pr 


o 




•k 














H- 










<1 




o 










rvj 




<1 










Z 










1 — 








sr 


VI 




w 




UJ 










H-l 






VI 




X 




X 




C 7 


X 




Q. 










VI 










X 


X 


CL 








>r 


3 




O 




X 






z 




X 






eH 


X 






X 




' 4 - 


3 








o 






a 










o 


* 


X 




o 




« 


< 




Z 




X 






H - 1 




< 






z 


• 


o 




X 






-> 




l-H 










»- 




X 






X 


rsi 


< 5 . 








1 










CL 






HH 










3 


'W 


X 




z 




W 


3 




O 




O 






(/> 




CL 


o 




o 


1 






a 




•R- 


3 




LU 




H- 






o 




a 


II 




X 




X 




H-l 




3 


< 




(X. 




O 






X 










cr 


-J 


o 




1 - 




-ft- 




o 


o 




<l 










X 






X 


*«r 






3 




-7 


V) 


►- 


h- 




X 






-5 




X 










z 




X 










CO 














o 








♦ 


X 




1-4 




X 


oe: 


o 




o 






a 




X 


O' 






X 


»-l 




X 




X 


u 


o 


UJ 




z 






o 




X 


• 




z 


X 


H 




1— 




X 


X 




CQ 




H-l 


• 




o 






f 




o 


—I 


<t 




V) 




♦ 










H- 












H- 




l-H 




X 




►— » 






z 


in 


>• 




X 


H 




X 




O 


a 




H- 




3 




a 




• 


a 


1 


< 




o 


X 


X 


X 




II 


• 




< 


1 


X 








•H 


•—4 


X 


X 








< 


h- 


Csl 


*-• II “5 


— • 




h- 


X 


o 




z 




•f 


h- 


• 






UJ 




ox 




•k 


-5“^ 


(M 




O 


o 


X 




< 






<I 


3 


4 , 1 ^ 




X 


— • • 




Z 


o 


•-<-5 






CL 


X 






X 


O 


-7 


3 


• 








•u> 


oo 


H-i 


MJ 


CL 








CL 


X 




rsi 


'U 


-R- 


3 


1 — 


X 




UJ 


•-4 O' 


oo 




•• 


CL 


rsjkO 




X 




X 




1 - 


•• 


• 


X 


3 






X 


• 


in>o 


X 


CM 


XXX 


V. •• 


W 


X 


II 


h- 




3 




CM 


O 


• 


u; 




h- 




^ •» 




II 


t— lUJLi^ 






1 — 






• 


LJ 


II 




X 










QCUJ 


-VO 


X 


^X 


OX 


II 


^o 




-J 


X 


o 


GQ 


3 






3 


X 




UJ 


w •rjww 


h- 


II 


XX^ 


1 X 


-h II 


X 




IN 4 


II 






It 


u 




h» 




Z 


\ XX 




ii 1 ^ 


ZJ-Z 




H- 


“7 


f-H 




< 


CN.H — 


z 


3 


u 


• 






X 


X 




1 


“7—1 




C# 


—1 


•wj 






1 3 




rvi 


0 "t 


C 9 


X 


II 




aw 




' X II ' 


—1 


X 


< 1 . 




X 


U 3 W 


X 




X 


II 


X 


-j II UU M XX 


p— 




X UJ 


xa n 


“7 V 


— j 


X 


f— • 


”7 




X 


II *7 


JC. 


u. 


UJ 


PsI 


LU 


• 7 X»-»«-<JSX 


CJL 


o 


XXX 




O 


Urn 


e- 


fNJ 


3 


3 f>l 


3 


f — 4 


z 




>— 










r-xo 






—1 




1-4 




V) 






VI 








X 






r- 




x<^x 






-M. 




z 


Ul 


VI 






Ul 




•J 




CJ 






L^l 




ox^ 






o 




»-H 




< 1 . 










•fr 




* 
















*fr 








-ft- 






-ft- 




♦ 




■Jl- 
















♦ 




•fr 




ir 






•ft* 






<w)OOIw) 


1^00 ooooo 


ooo ooo ooo 


ooo 



57 



C**C^LCU/iTE THE PARTITICN FUNCTICN I 

^ 2 2 = 2-»ZJ( L »*GJ(U CHI00580 

0*EEGU THE CALCULATICNS .1 1/CH FRCH THE FIRST FBECUENCV 



oooo 

U> 0 *-«CNJ 

IOnOMJnU 

oooo 

oooo 




oooo 



oooooooooooooooooo 

fT) in r* ou 0^0 n) nt u T >o 00 of>o r-H 

'i->r*r-r^f^f^n-r*r*-Goaooo 

oooooooooooooooooo 

oooooooooooooooooo 

>■^1 M ' M n <1 M-itI < 1 < ! ■<>— n 



OOOOOOOOOOOOOOOOOO 



o 


o 


ooo 


I'l 


nT 




cu 


UU 


« 3 COf-< 


O 


o 


oo^ 


O 


o 


ooo 




1-4 


1 >1 II i 


X 


X 


XXX 


o 


o 


ooo 















CM 












sr 






w 












O 






CO 












'4- 






GC 


















< 












• 






O 












»-♦ 


















1 






— • 






























* 


00 




m 












-> 






W 














o 


f— 


cn 
















►-4 


X 




tM 








CL 




X 


< 














O 




o 




1-4 








ill 


O 


o 


+ 




3 








1 




•-4 














• 




—1 


^4 




— * 










— • 


o 
















CM 


UJ 


VI 




tM 












<c 


QC 














1/) 


OL 


< 




CO 








1 ^ 


QC 


< 


O 




X 








I — 


< 


GL 






< 








•UIX 


o 




V. 




o 








-400 


• 


h- 


CM 




+ 








1 cr»2 


UJ 


2 


* 












^ •LU 


-J 


H-l 


* 




f— 4 










• 


O 














-*Hr ^ 




X 




>• 


VI 










CO 




Ol 


w 


X 








<1—1-5 




LU 




2 


<1 










CO 


X 


00 


X 


o 








<-30 


QC 


OL 


QC 


3 










ooo 


< 


X 


<L 


o 










HJJ* 


o 


1— 


O 


LU 


—2 








3QCX 


• 




1 


QC 










*«-xa 


a 


< 




X 


fOX 








3*2 


2 








-^< 








1 — Ol 


< 


>• 




o 


LOX 










• 


CO 


CO 


2 


QC3 








•fr ww 


— % 




CL 


•-4 


< •. 








* -3-1 


p-4 


UJ 


< 


CJ 


o— 








— 0<A 


w 


o 


o 


2 


1 2 








— I*UJ 


CO 


— J 


w 


O 






»-4 




w«^a: 


QC 


< 


* 


X 


f-4X 




X 




-3-Jw 


< 


> 




CO 


^^< 




o 




O-'l- 


O 




• 


UJ 


COX 








111 -5 a: 


o — • 




GO 


OL 


Q^»-l 




o 




cLisiar 


20 (Mcn H 


< 




OL 


<x 




z 


o 


x*co 


hh2 ww o 


LU 


• 


a 


oo 


• 


1-4 


o 


0'-'O*»'-3 


3*“* I/3CO • 


X 


in 


o 






I— 


o 


sOXi-irv|»-t 


*3 CXQC»-<— 








*o 




<A 


o 


••-ICO* X 


• ^ <<XOJ 


lU 


II 


UJ 


IX'IO 


-^ooo 


o 


•» 


^GlOI— o 


CM3 OOO^ 


X 




X 


•CM O 


It tl II 


o 




tl A* ^ 


HI ^ II 11 H VI 


1 — 




>- 


tl vnuo 




o 


11 z 


-IOCL>vi-i03 II 3 — — (Xm 








-^NOorn 


o-^cvjrn 


-J 


5ti-i 


• II f\jx2 II — II •HCMrn< 


LU 


X 


UJ 


2— '2 


UJwww 


< 


JIO'TX II *0 




2 


< 


2 


'-'LU *-<0 


tX.COl/JVI 


o 


iM ! 


n 0 * 5 * II 1— *-i^nicoi/3co^^--r 




X 


►-4 


XI— f-'t— 


u-oiQCa: 








X 


►-4 


X 


<•-•2 CJ 


tl 




LJ It 




X 


X 


X 


XXLJLJX 


j&WWW 




LJ JtLJ LJLJH-OLJ Jfc XOLJLJ O^ 


UJ 


o 


UJ 


330 CJ 3 X 




o 




f-l 








h- 




r- 


ULi 






UJ 


UJ 




UJ 


tM 



o 



•Ifr 

ooo 



LJ UJ 

* * 

♦ * 

ooo ooo 



58 



FPOGRAM TO CALCULATE THE UNCERTAINTY IN 
























Z 


(X z 




















3 


X— 3 










































3 


3Z3 






• 














w 


Z^3^ 






UJ 














> 




UJ 




LJ 














3 




o 


ut 
















Z 


1 — >^x 




IX 
















X 


33 z 


<l 


3 


CJ 














Z3 


<izxx 


J. 


H- 


CO z 




1 




CO 






33 


crxocr 


LJ 




h- U 




1 




1-^ 






H-lX 


X3X3 




□C 


z z 




1 






z 




H-C3L-. 


>-<3a:z 


UJ 


UJ 


< 




1 




Z 


— < 




oxz 


ZXXh-i 


X 


a. 


h- z 




1 




0 


z a: 




X (3 


i-ia: 


H- 


X 


CO 0 




1 




1-^ 


u e> 




C03V. 


XCO> 




UJ 


Z CD 




1 




H 


V. a 






Z X3 




h- 


OQC.0C 




— ^ 1 




<1 


Z *-• QC. 




COZ^ 


33ZZ 


i-i 




LJLJ<I 


t 


0 1 






I-I W> UL 




C03h- 


•"Kll— IX 




CJ 


XLJ 




^ 1 




3 


> 




3»— iZ 


H*Zh“3 


LJ 


r- 


> 


C3 


Ul 1 




W 


3 Z UJ 


X 


QLh- -a. 


<IO 3 


UJ 




OCh-iLJ 




0 1 




3 


Ul C3 X 


0 




OO 




<xz 




a I 






X n- 




C^COXXKZQC 


<l 


0 


<00<J. 


o: 


C3 1 


»-«cMrn<r 


C3 


COUI 


CO 


ZXZ3I-I<1<IX 


ou 


Ul 


CO 


Ul 


1 






couico CO 




<II--I03>cxh- 




a. 


Ul 2T 


•-4 


1 




< 


ujQCcr CL 


CO 


X>33 1 XCOX 


K 


CO 


LJUIUI 


u. 


1 






UIUIUI-<»0 


•-4 


<C 


ZZH-IZ3 


Z 


UJ 


UIH-O 




1 




CC 


a:x>e?h- 


CO 


0 : 0:3 


o>oz 


UJ 


oc 


z<o 


n- 


1 




0 


OXZ'^CO 


z 


<<X*-i 3X 


X 




301 


z 






X 


UlCO*-» 


0 


33Z3H-3 QC 


Ul 


X 


UI3C-J 


Ul 




Ul 




30 Ul • 


0 


<l 3 0 <t <I ►-•O X 


oc. 


H- 


xo> 


Q 




0 


K 


zzao 




h-0*-»Xh- Z^X 


3 


►H 


I-3X 


»-4 






3 


zn-i-io 


<t 


oxn- 


03Z»-i 


CO 


JC 


<C •* 






X 


X 


C3X 


h- 


1 — — l<aU.(X2:QO 


<l 




00 <3 Z 


z 


1 


z a 


Z 


0 




ot-u OCX 


Ul 


Q 


<2 Ul 


< 


1 


Ul z 




-*z<o: 


Q 


xxo 


X<QCX 


3 : 


UJ 


XO<3 




1 


zozo 




1 — ►iZUIUI 




X 0CZXX<X 




q: 


h-O 


CO 


1 UIOUIZ 


Q 


^ <H1Q^ 


3 


h- X 1 


H-l|— XI- 


UJ 


< 


Z a: 


< 


CO 1 


LDOHO 


Ul 


X— 00-3 


<t 


XXX -tz 


q: 


3 


on- 


X 


<t 1 


>»-az 


<r 




Z 


t«|^XcO*«X< •• 


3 


<3 


er 3»-4 




CJJ 1 


x*-icrcj 


H>4 


3^xn-< 


OC 


0 1— 


XO 3 


f- 


CO 


e>3Z 


CO 


1 


C3ZCJCQ 


3 


H-uji— zer 


X 


•• 


••XX ••z 


<i 




o<t •• 


< 


1 


>cr 


0 


x^QCOUlUJ 


I— 


CO>**l 


Mi3<XX»-i 


QC 


>—• 


cezz 


e? 


1 


x< 


Ul 


CT3»-*3X 


z 


Q3xcr<xx3 


UJ 


X 


XQCUJ 




1 


V3 


OC 


UJCOjCH-12. 








Q. 


u 


uie> 


X 


1 






Q.C0UJ Ul 








2. 




uiH-> 


u 


1 




Ul 


ZUIZCOI— 


X 






UJ 


z 


XZX 


< 


1 




X 


Ul QCH-K 


X 






J— 


1-^ 


^-*-•13 


Ul 


1 




h- 


1— X3C3<1 


h- 







I 



LJ O W O O O O e? O O O O CJ O O O U? o CJ o ^ U? CJ o o e? O CJ 



59 



ooooooooooooooooooooooo 
^ csj m u su CO o rvj o I ^ m vu CO o ^ ^ * 

O OO O OO O O •-< iM ^ ^ rH f-4 1-4 IM INJ <\J l\i 

00000000^00000000000000 
OOOOOOOOOOOOOOCPOOOOOOOO 

M M n iHHhhI II < 1 M H * 

oooooouuooooooouooooooo 



0 


•» 


0 


^ • 


0 


0 -< 


nO 


0 ih 


f -4 


• ST 


•» 


•• m 


IM 


m •• 


w 


m 


:> 


0 0 


<i 


• 0 


CO 


CM •• • 




^ CA ^ 


X 


- OCT 


0 


fH -orn 


•» 


«u> 


— • 


mo »KM 


rH 


•om • 


nO 


0 «»Hm 


w 


(JO •»GO«^ 




moo • ** 


0 


t- 4 Goom 


•» 


•^CMCM 




cr •r-«U 9 




^ • 




mf- 4 cncp 


CO 


t-H ^m^u 


oc 




< 


•00<M 


0 


•• «o •* 


•» 


#cD^m 




om' 4 “»H 


m 


•.PI •MA 


w 


r-cMQur* 






X 


in^^o 


•* 


'Tcor- • 




^p-o •• 


0 


• •H • • 


>o 


*- 40 P)cr* 




•. • •*( 7 > 


—*-5 


m ^m •• 


00 


0 • 


OJOi 


UJ'H •.CM 


*- 4 Ct 


vU ••'TrH 


wU. 


mr->OPi 


X •* 


moo 00 


<x— 


•cr* • • 


XO 


•HCPO^H 


-S^O 


«. %\0 m» 


•-W 


sT *-4 


—-5 


cT) •-mo 


Of^ 


O^-OUJ 


CNJ •• 


0 •uico 




•om-r 


a-'O 


•-mm^H 


ZXt -4 


(Mm»H • 


OJ<L •• 


0 €04 #<M 


QXO' 


•H •»'T — 




cx^u' -m 


XXCL 


ujcLmmo 


UJOO 


ooooo 




uu 0 • • 


Q 


1 — <X • - 


X<1<1 


•.PUM 


LJUJUJ 




OGIUC' 


-• a^pi»H^ 





CO 

• 

m 


• 


fH 


- - OJ 


Oi 


— OC, 


•> 


- 0 


• 


xcoc/^-*co 


M 


OIOI- CO 


• 


COOiOiXOi 


OX 


ojodcx a: 


LU 


QJOOUOX 


xco 


crujmuj 


XUJ 


OCJUIUI 


^cu 


uu a CO 


VJCX. 



>ru) Oi. o 

• • • LUO OOi 

o •• ♦‘Oh- ZQ 
oo < <- 

LL. •. t • - •— » 

••Otntn* ^ Q»H» 
« •LLLL *0 OI 
Qin •* -o • Q^oo 
uju.» • «in I 

•• mu- ou.H -1 
O- l-h-LL •• o ••X 
<<t — c/)- O 

II Q 



ro 

o 



ex 

COCOO UJUJ o 
• ••-*2:001— U- 
— mir\ 3 aioo< 
csji^iF-iaiooia: ^ 

*-iujaiz>Ho -»xo 

•• •* •»«-hX*^>* 0 - '^O*— *0 
orom_io 2 xo "vQcoi— 

• I M «» - - - • - - 

U-XXX^X^XXXXXX 

rnrH*- 4 r-*»— 



a z 
I <•-• 
i-iXU- 







UU UL tJU kL LJL ex. a. UL UL CL UL OC. 
OOOOOOOCJOOOOO 



WOOOLJOOOOOOOO 
00 00 000^0 

»HC\j^m u ' m >ur-co CO cn 



0000000000 
>r m sor- coo o*-«osjn> 
cvjcNj iMOj cgrvj m 0 1 n 1 p I 

0000000000 

0000000000 

<1 <1 M l» II |l l» .IHH 

o 5 cj 5 oo 5 o 5 o 



00 

Pirn 

00 

00 

* M "H 

XX 

00 



0000 





























0 




















X 




















q: 










1 - 










< 










3 L 










OC 










0 










1 - 










0 










< 










X 










0 










Z 










UH 






h-l— 1 — I— 




1 - 










X 




</) 


0 0000 




0 




LU 


- ^CMPInT 








-J 


<i mmu>m 




X 




OJ 


y •»•»•»•» 








< 


c vO'O'OO 








•-H 


CD wwww 


Q. 


z 




Qd 


•. LULU LULU 


0 






< 


a- h-Hi-i- 


h- 


X 




> 


•. t <1 


t/> 


CO 






1 — ^QCiQCcxa: 


CM 




•H 


1 — 


— Ojc JC-S 3 : 




X 


'U 


- Z> 


00-*-^— ‘ 


•H • 


X 




a. 


or^^HCMm^ 


MO 


1 - 


*-4 


z 


•-4 •. • • • • 


^ • 




II • 


l»H 


-vOOc^cjrcDf 


xo 


X 


^f -4 




mwoiaiujuj 


UJ 


fvi 


II 


LU 


wUJ • « • • 


m • 




> 0 --* 


X 


UJ 1-0000 






wl 


1 — 






0^ 




UUCXCi-CLU-U. II 


II CJCi. 




o-> 


LJ 


IX ^ *-i *-1 *-H 2 . Cj 


h- 


0 


<L 










UJ 


0 




z 


'U 


' UL 


0 








* 


m 








■th 






*#«• 




000 


000 



o 



60 
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C♦*^SSU^E A 6CLT2MAN DISTRIBUTION FOR EACH LEVEL 
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C**I^CPE^ENT THE TEMPERATURE EY ONE DEGREE THEN CALCULATE 
C ALL THE VALLES OF CHI AGAIN 
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